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Systems of nonlinear ordinary differential equations are constructed, for which 
the general solution is algebraically expressed in terms of a finite number of par- 
ticular solutions. Expressions of that type are called the nonlinear superposition 
formulas. These systems are connected with local Lie groups transformations on 
their homogeneous spaces. In the presented work the nonlinear superposition for- 
mulas are constructed for the case of the SO(3, 2) group and some aspects in the 
general case of SO(n + 1, n) are studied. 

1. Introduction 

> . 

Q\ . Let us consider a system of n first-order differential equations 

g: d?(t) = X ^(x\x 2 } ...,x n } t), (i = l,... ,n, (1) 

o 

where the dot denotes differentiation x M (t) with respect to time t. It is known for 
a very long time that, in some cases which we will specify later, it is possible to 
express the general solution as a nonlinear function of a finite number of particular 
£H \ solutions; it is of the form 



F(xi,x 2 , . . . ,x m ,ci, . . . ,c s ) , xGR", (2) 



■ where xi, . . . , x m are particular solutions ([[]), ci, c 2 , . . . , c s are arbitrary constants, 
and x(t) is the general solution. These relations are called nonlinear superposition 
formulas. Here x means a vector with elements 

An example of such systems is a homogeneous system of the first order linear 
differential equations, in which the general solution is expressed as a linear combi- 
nation of n linearly independent particular solutions. The other known example is 
the Riccati equation 

± = a(t) + b(t)x + c(t)x 2 , (3) 

where a(t), b(t), c(t) are continuous differentiable functions with respect to t. In 
this case, for any four solutions Xi(t), i = 1, . . . , 4, the relation 

xi{t) - x 3 (t) x 2 (t) - x 4 (t) _ Ui - u 3 u 2 - u 4 ^ 
Xi(t) — x±(f) x 2 (t)—x 3 (t) U1—U4 u 2 - u 3 ' 
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where Xi(0) = Ui are initial conditions, is valid. 

In the general case, such systems of differential equations are connected with the 
local Lie group G of transformations on a factor space M = G/Gq, where Go is a 
Lie subalgebra of G [ffj. We recall this connection briefly. 

By the local Lie group G of transformations on the M, we understand a smooth 
mapping ip : G x M — > M, (we use the abbreviation (p(g,u) — g ■ u), for which: a) 
e ■ u = u, for any u G M, where e is the unit element of the group G, b) for any two 
elements g±, g 2 G G and any u G M is g 2 ■ (<?i ■ u) = (<72<?i) " M an d c) g ■ u = u for 
any u G M imply ^ = e [0. 

In the local coordinate system, we write x = g ■ u as 



x 



»-r(a\...,a N y,...,u n ), M = l,...,n, (5) 



where is a dimension of the group G and a r , r = 1, . . . , N, are their local coordi- 
nates. For x M (a, u), we can write 

|: = Ee(xK(a), (6) 

n ^ 

where the vector fields XJx.) = — ^'(x)^; — obey the equations 

~, ox^ 

[x r ,x s ] = f;4x t , 

in which structure constants of the Lie algebra of the group G. 

Conversely, the vector fields Xfc(x) determine the infinitesimal action of the local 
Lie group on the space M uniquely [§]. 

Let g(t) = (a 1 (t), . . . , a N (t)), t G R, be a curve in the Lie group G such that 
g(0) = e. This gives a curve in the space M. Derivation of equation (|5|) with respect 
to the parameter t gives, by using (^]), the system of differential equations 

^ = E^(x)Z r (t), /i = l,...,n. (7) 

r=l 

In our paper we will deal with such systems connected with the Lie group SO(n + 
l,n). 

If the system of equations (|I|) has form ([/]), then there is a curve in some local 
coordinates on the Lie group G, which acts on the factor space M. In this case, it is 
possible to find the superposition formula |TJ . Any particular solution of the system 
(0) can be written in the form 

x fc (t) =g(t)-u k , (8) 
where = x^(0) is the initial condition. 



2 



Xl 


= f(a\.. 


.,a N ,u 1 ) 


X-2 


= f(a 1 ,.. 


.,a N ,u 2 ) 




= f(o 1 ,.. 


■ , a N , Ur) 



We express the action of the local group G by using action of this group to few 
points of the space M, which is supposed known. In principle, it means to find from 
the system of the equations 



(9) 



the coordinates of the group a 1 . To be able to find the group coordinates a 1 , i = 
1, . . . , N, we should use the action of r points. It is evident that the number r must 
fulfil the inequality nr > N, where n is the dimension of M, and N is the dimension 
of the group G. 

If we solve this, then we are able to express the elements of Lie group G by 
means of the known transformations of r points in the form 

g = #(Xi, X 2 , . . . , X r , Ui, u 2 , . . . , Ur) , (10) 

therefore 

x = 5f(xi,...,x r ,u 1 ,...,u r ) -u (11) 



holds. Now we see that formula (|11|) is invariant with respect to the action of the 
local group G. 

If Xj(t), i = 1, . . . ,r, are known solutions of the system (|7D for given functions 
Z r (t), then any other solution x(£) of that system is given by (|TT|) . Therefore the 
relation (|TT| ) is the superposition formula 

For example, the Riccati equation (Q) is connected with the Lie group G = SL{2) 
that acts on the space M as 

/,x ,,x a 2 x(t) + ua 22 (t) 

git) ■ u = x[t) = — r — , 

W W a u (t)+ua 12 (t) 

where we represented the elements of the group G by matrix I 041 fll2 ] , with the 

\ a 21 «22 / 

determinant equals to 1. 

In the eighties, there where a lot of papers which were devoted to the 

systems of differential equations of this type and especially to finding the superposi- 
tion formulas for these systems. For the most part, the authors study the systems, 
which are connected with Lie groups SL(n, R) or SL(n, C). In the general paper [] j 
the methods are described that are used for constructing the superposition formu- 
las. They are demonstrated with examples, in which the Lie groups SL(n, R) and 
0(p+l,n — p + I) act on simple projective spaces. 

In P] the more general case of the projective matrix Riccati equation are studied 
for SL(2n, R) and SP(2n, R). In paper Q, the systems that arise from the action 
of SL(n, C) on the factor spaces SL(n, C)/0(n, C) and SL(2n, C)/SP(2n, C) are 
studied. Paper ||, is devoted to the superposition formulas for rectangular-matrix 
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Riccati equations on the space SL(n + k, C)/P(k), where P(k) are special maximal 
parabolic subgroups of the SL(n + k,C). In ||, the authors deal with systems 
connected with the Lie group SU(n,n), and in 0, th same method for SO(n, n) 
is used. Finally, in paper ||, the authors study the systems of equations that are 
connected with the action of Lie group SL(n, C) on the space M = SL(n, C)/Go, 
where G is a special non-maximal parabolic subgroup. 

The authors mostly used a set of special solutions for reconstructing the group 
action on the space M. This approach simplifies the solution of @. But on the 
other hand, we are not able to use the resulting superposition formulas directly for 
constructing the solution of the system on the basis any set of particular solutions. 
By means of the action of groups elements, we must first transform our particular 
solutions to a special set of particular solutions used in superposition formulas. 

In the next, we study the systems of equations (0) which are connected with 
the action of Lie group SO(n + 1, n) on the space M = SO(n + l,n)/P, where P 
is one of the maximal parabolic subgroups. To our knowledge systems of this type 
were not studied so far. In distinction to the papers cited above, we do not choose 
a special set of particular solutions for reconstructing the group action on M, or, in 
the terminology of paper , we will construct group invariants which give nonlinear 
superposition formulas in implicit form. 

2. Lie group SO(n + l,n) and its Lie algebra 



In this section, we fix the notation. The Lie group SO(n + 1, n) is a group of real 
matrices G with dimension (2n + 1) x (2n + 1) that fulfil the equations 



G T • a ■ G = a , where (7=001, (12) 




G T denotes a transposed matrix and I is a unit matrix with dimension n x n. 
Matrix G is written in the form 



511 


T 


T 
gl3 


g21 


G22 


G23 


g31 


G32 


G33 



G = g 2 i G 22 G 23 , (13) 



where g denotes a column vector, and g T is its transposed vector represented as 
row. Therefore gg T is an n x n matrix, and g r g is an inner product. 

If we insert (13) into ([12]) we obtain the following formulas for elements of the 
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matrix G 



9u + 621631 



531S21 



1 



ffllgl2 + G^gsi + G^ 2 g21 = 
01lgl3 + G^g 3 l + G^g2l = 



>12gl2 + G22G32 + G3 2 G 2 2 — 



(14) 



gl3g?3 + G23G33 + G33G23 







>12gi3 + G2 2 G 33 



GX,Go3 — I 



r 32 VJ 2 3 

In this realization the Lie algebra so(n + 1, n) is given by real matrices 

/ x T z T 
A= -z H W 
V -x Y — H T 

where W T = -W a Y T = -Y. 

In the group SO(n + l,n), we take a subgroup Go that is generated by the 
matrices 

1 
D 
(D T ) 



G = 



where the equality Z + Z T + zz T 
represented by the matrices 



/ 1 z T 
-z I Z 

V i 

holds. The factor space M = G/G can be 




where X + X T + xx T = 0. As coordinates on M, we choose x and the antisymmetric 
part of the matrix X; this means that 



Y = X + - xx 



(15) 



The action of the group G on the factor space M in these coordinates can be obtained 
from the equation 

\ 



T 
gl2 



/ 9u 

g21 G22 
\ g31 G32 




/ 1 



V 




1 





z 

Z 

I 



If we compare the coefficients on both sides of this equation, we obtain the following 
formulas 

D = G 22 + g 2 iu T + G 23 U 

x=(D r )- 1 (g 1 2 + ^ 11 u + U T g 13 ) (16) 
X = (G 32 + g 3 iu T + G33UP- 1 
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for the group elements. Because we restrict ourselves to the local Lie group, we can 
suppose that matrix D is invertible. 

Starting with the action of this group on the space M and by using the expansion 
to the first order, we derive an explicit expression for vector fields, in the basis of 
the algebra so(n + l,n) in the representation (T g f)(m) = /(g^ 1 ■ m). Specifically 
we get 

d 



Y 
1 u 



dY l3 



B 1 ™ B 
Xi i ^ -2^2. 



dxi 2 f^ x dY ri 



d - d 



2 I Br 2 ^ BY 

r= l \ ^ / L/^ r i r s=l ul rt 



n B n 8 



r=l vX T r )8 —i OY rs 

where we define Y rs = —Y sr . 

The system of differential equations (0) is in this case of the form 

x(t) = a(t) - B(t)x - Yc(t) + -xc T (t)x - YC(t)x 

Y(t) = A(t) + i(xa T (t)-a(t)x T )-(YB(t) + B T (t)Y)+ (17) 

+ ~ (Yc(t)x r + xc T (t)Y) - YC(t)Y 
2 

where a(t) and c(t) are vector functions differentiable with respect to t; A(t), B(£) 
and C(i) are differentiable matrix functions, for which A T = —A, C T = — C and 
Y T = -Y. 

This system of differential equation is the same as in the paper ||, where is given 
classification of all systems of nonlinear ordinary differential equations of this type. 

3. Representations of the action of the group by means of particular 

solutions 



The system of equations ([171) arises from the action of the Lie group SO(n+l, n) 
on the factor space M = SO(n + l, n)/G . Therefore, one can find the superposition 
formula. The action of the group SO(n + 1, n) on the space M is given by the 
relations ([H]). We try to express this action by means of any known solutions 
X*. = g ■ Ufc. If Xfc(t) are solutions of the differential equations fllTf ) with the initial 
conditions Xfc(O) = u^, we obtain the general solution x(i) with the initial condition 
x(0) = u in the form 

x(t) = G(xi(t), . . . ,x r (£),ui, . . . ,u r ) • u, 
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where g(xi(t), . . . , x r (t), ui, . . . , u r ) = g(t) is the expression of an element of the 

group G in terms of the known transformation. 

Dimension of the group SO(n + 1, n) is equal to n(2n + 1) and the dimension of 

fi{fi -)- x) 

the space M = SO(n + l,n)/G is . So, to express iV coordinates of the 

group elements by using the known solutions, we must know at least r solutions, 

where r fulfils the inequality n(2n + 1) < r. Consequently, for n — 1, it is 

sufficient to known three particular solutions, and for n > 1, we must know at least, 
four particular solutions of the system (|T7D. 

By X, we denote a vector, and by Xj = Yj — — Xjxf , a matrix for and i th solution 

of the system of differential equations (|T7| ) with initial conditions Uj and U,; Dj is 
an invertible matrix defined for an appropriate solution by the first equation in ([T~ 



Further, we denote Ujfc = u L — u^, x^ = x, — x^. Similarly, we define Xjfc and XJik- 
Temporarily we suppose that all matrices which we will use in our calculations are 
invertible. 

From equation ([[(]) we obtain some coordinates of elements of the group G. 
Readily we discover that for indices i and k, the relations 



813 = 

G23 

G33 

gia ~- 

G22 

G32 



Dfx fe - 

.TVTT-l 



mr • (Dfx 

= (Di - D fc - g 2 iu| fc j 
= (X,D, - X fc D fc - g 31 uf fc ) • U£ 
Dfx, - g n ui - Uf gi 3 
= D, - g 21 uf - G 23 U, 
= X 4 D, - g 31 uf - G 3 3U, 



are valid. If we now apply these equations to Q14|) , after a simple algebra we obtain 
the following formulas 



9u + giifei + 831621 



1 



Df ■ (g n xi + g 31 + Xf g 2 i) = ^ 



D 



x,xf + Xf + X fc ) • D fc = muT + Uf + U fc 



(18) 



From the second equation in ([Tj|) we can express 



g21 
g31 



-(xf fc )" 

Xf (X^.) 



(DfrV 
- fDfr 1 



fx 



T\-l 



U; 



Xf (Xfcj 



S'nXfe 
1 



(Df) 



-1 



flfiiXfc - (Df ) u fc 



(19) 



Next we use the notation 



^ik 

hi = gnXi 



Xjxf + Xf + X 



Ui uf + Uf + U 



fD 



U; 



(20) 
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It is easy to see that the relations 



XjX^ + Q ik + Q ki = and 
UjUj[ + Ui k + io k i = and 



T 



a; 



are valid. By simple algebraic calculations from equations 
j, and k, we derive the following formulas 



and (|l^) for any 



Uik 



V = 



<?1 2 1 + 



x^r^+x.Tfxi.r'h 



-hfXr^X 



L h, - h^x^^xf,)- 1 ^ = 1 



(21) 
(22) 

(23) 



We found the system of equations from which it is possible to determine the matrix 
(|T3|). Since we restricted ourselves to the neighbourhood of the point t — 0, we 
can suppose that matrices Dj(t) are invertible. From the relation (|21~1) it follows 
that in this neighbourhood matrices flik{k) and u>ki are simultaneously invertible or 
non-invertible. So, it is easy to see that, if matrix uj^ is invertible, it is enough 
to know only one solution of this system Dj for one value of i. We can find other 
Dfc from equation (^l|). If we know D« explicitly, the equation (|2l 



is a quadratic 

equation with respect to g\\. The odd equations can then give the matrices Dj. 
We enunciate the previous account in terms of the following theorem. 



Theorem: Let Xj(t) and Yj(t), i = 1, 2, 3, be three solutions of the equation 
(|T7|), Xj(t) = Yi(t) — -Xjxf, Ui = Xj(0), Uj = Xj(0) and let the matrices XJik 

be invertible. Then, there is a neighbourhood of the point t = 0, matrices Dj(t), 
i = 1, 2, 3, Dj(0) = I, and function gn(t), gn(0) = 1 such, that, for them formulas 
(H), (II), and (|D are true. 



4. The superposition formulas 



In the previous sections, we formulated conditions (pTj) (P3|) that are valid for 
any three solutions of the system (|TTD. Although we are not able to solve the system 
(|2l|)-(|23D explicitly, it is possible to derive, from them, certain relations for solutions 
of the system of differential equations. 

We suppose that we have five solutions of this system, and there exists the inverse 
matrix to io^ for any i ^ k, i, k = 1, . . . , 5. Then from the equation (pl|), we obtain 

If we now multiply the equation ( pT| ) by this equation from the left, and then, by 



(|2l| ) from the left for the couple (rz), we obtain equation 

(24) 
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that is true for any i, k, and r. We multiply this equation further from the right by 
the inverse equation for the ternary (sti). Then, we obtain 



which is true for any (ikrst). If now we put s = k in this equation, we obtain, for 
any four solutions of (|ITD, the relation 



^-QjQskQjQri&i = UJ^UskU^LUri . (25) 

It means that matrices 0,"^ fl sk Q~ k Q ri and lu^ u a yjj~ k uj r i are similar. Therefore, all 
their invariants are identical. 

Now, we will use this interesting property of the solution of differential equations 
(|T7| ) to obtain superposition formulas of the system (|17D for small n. 



For n = 1 we get the Lie group SO(2, 1). In this case, the vector x is reduced to 
the number x, and Y = 0. From this we have Xj = — x\ and flu- = — (xi — Xk) 2 ■ 
In this case, equation (p5|) has after the extraction form 

tlC g t> 30 y 00 Ui g %ii Hi •<{• Hi £ 



dC g "JO 2 *jC y iC tjC fj^ Hj 2 r Lii f Hi 

That represents the well-known superposition formula @) for the Riccati equation 
(0). This is a consequence of the isomorphism between so(2, 1) and sl(2). 

Now, we will study the case n = 2, that is the Lie group SO (3, 2). In this case, 
we have 

Y=( ° x *) and X=if ~ X \ ~^ + 2x 3 \ 
\ -x 3 J 2 \ -xix 2 - 2x 3 —X2 J 

If we denote the components of the i th solution by obtain 

det(U lfc ) = \ (4(4 l) - uff - {ufu? - ufu^f) 
and the determinants of the matrices and uj^ are 

det(Q ik ) = ^A 2 ik and det(u; ifc ) = - 8% , (26) 

where 

A o (*) o ( k ) i (i) (fe) W (fc) i <r (i) n (fe) . (i) (fe) (i) (fc) 

Aj fc = 2x3 — 2x 3 + x 2 x\ — x\ x\ and 0^ = 2u\ ' — 2u\ ' + u 2 u\ — u\ u\ . 

(27) 

We see that the conditions for the matrices \Jik to be invertible imply the invertibility 
for matrices uj^- 
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If we now take the determinant in equation (|25|) we obtain for any four different 
solutions the equality 

A s fc A r j 5 s k $ri /<->o\ 
A - ' A = "7 ' 7 ' \ 2H > 

t-^si ^rk "si "rk 

As it was mentioned above, in this case, we can construct the general solution of the 
system (|17|) by using four particular solutions. Take now five solutions, for which 
from (|28|) we obtain independent equations 







A 24 Ai 3 


_ #24 


#13 


A 34 A 12 


_ #34 


#12 






Au A 23 




#23 ' 


A 14 A 23 




#23 


A 25 


A13 


_ #25 <>13 


A 35 


Ai 2 


_ S 35 5 12 


A 45 


A 


A15 


A 23 


^15 #23 ' 


Ai 5 


A 23 


5l5 #23 ' 


A15 


A 



(29) 

12 _ #45 £l2 
\2A ^15 ^24 

Equations (SO) are understood as a system of linear equations for o , and 

x^ . This system has a solution when its determinant D(t) is different from zero. 
By direct calculation for t = 0, we obtain 

.D(O) = (5l5-525)(5l3-^14 + 534)-(^15-^35)(5l2-5l4 + ^24) + (5l5-^45)('5l2-5l3 + ^23) • 

(31) 

(i) 

The terms in the second parentheses do not depend on u 3 and are not identical zero. 
Because the terms in first parentheses depend on —u^ \ this determinant is not 
identical zero for any possible initial conditions u, and Uj. As the determinant D[t) 
is given by the solutions of system ( |TTD which are continuous and -D(O) ^ 0, there 
is any neighbourhood of t — 0, in which the determinant D{t) ^ 0. We see that, 
in this neighbourhood we can determine, from the system of the equations fl3"0]). 
the solutions xf\ x 2 an d of the system of differential equations ( |TTD by using 
the particular solution Xj, Yj for i = 1, . . . ,4. In other words, formulas (p8|) give 
the implicit form of nonlinear superposition formulas for the system of differential 
equations (|17D which is connected with the action of the Lie group SO(3, 2) on space 
M. 



Comments: Equations (^) imply that four solutions are not fully independent. For 
example, if we know three solutions Xj, Yj, i = 1, 2, 3, and from the fourth we know 
x^\ we can obtain and x^ from fl2"9"|). This is a consequence of the fact that 
the reconstruction of action of the group requires only 10 independent functions. 

5. Conclusions 

The main results of this paper are the following: 

1. We constructed systems of first-order ordinary differential equations that 
arise from the infinitesimal action of the local Lie group SO(n + 1, n) on the factor- 
space M = SO(n + l,n)/P, where P is one of the maximal parabolic subgroups of 
SO(n + l,n). These systems allow a superposition formula. 

2. We found any set of invariants for these systems which are expressed in terms 
the solutions. These are invariants of matrices VL^. The matrices VL^, play, in our 
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similar role as the matrix anharmonic ratio for projective matrix Riccati 
equations in ||. 

3. In the case of S0(3, 2), we proved that, from these invariant it is possible 
to find the general solution of our system on the basis of four particular solutions. 
Therefore, in this case, these set of invariants gives implicit nonlinear superposition 
formula. It is necessary to note, that, even though the local groups SO(3, 2) and 
SP(4, R) are isomorphic, our system of differential equations differs from the one 
studied in for the Lie group SP(4, R), because we use another maximal parabolic 
subgroup for constructing the space M. 
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